Quantum computing in optical microtraps 
based on the motional states of neutral atoms 
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We investigate quantum computation with neutral atoms in optical microtraps where the qubit 
is implemented in the motional states of the atoms, i.e., in the two lowest vibrational states of 
each trap. The quantum gate operation is performed by adiabatically approaching two traps and 
allowing tunneling and cold collisions to take place. We demonstrate the capability of this scheme 
to realize a square-root of swap gate, and address the problem of double occupation and excitation 
to other unwanted states. We expand the two-particle wavefunction in an orthonormal basis and 
analyze quantum correlations throughout the whole gate process. Fidelity of the gate operation is 
evaluated as a function of the degree of adiabaticity in moving the traps. Simulations are based 
on rubidium atoms in state-of-the-art optical microtraps with quantum gate realizations in the few 
tens of milliseconds duration range. 

PACS numbers; 03.67.Lx, 32.80.Pj, 42.50.Vk 
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I. INTRODUCTION 

The development of tools to prepare, manipulate and 
measure the quantum state of a physical system rep- 
resents one of the great challenges of modern science 
and, in particular, it is essential for applications in quan- 
tum information processing such as quantum comput- 
ing. At present a few systems have been identified that 
should permit quantum computation: molecules in the 
context of NMR R] , ion-traps |^ , cavity QED with pho- 
tons and atoms H, solid state devices such as quantum 
dots 0, ||, H , and trapped neutral atoms ^, ^, |l^, ^ . 
For Rydberg atoms in high Q cavities the engineering 
of entangled states and the implementation of quantum 
logic have been demonstrated [Q, furthermore a quan- 
tum gate has been performed between the internal and 
external degrees of freedom of an ion in a trap jl^ . In 
NMR systems, quantum algorithms on a few qubits, e.g., 
Shor's factoring algorithm, have been reported pH |. 

Neutral atoms are promising candidates for quantum 
computing for at least two reasons: (i) techniques of cool- 
ing and trapping atoms are by now very well established 
|l5[ ; and (ii) they are comparatively less sensitive to deco- 
herence, e.g., interaction with the 'classical' environment. 
Neutral atoms can be stored and manipulated in optical 
lattices J|l6[, standard dipole traps and microtraps 
p8| , p^ , pO[ In particular, magnetic ||l8| and optical 
microtraps ||19| , ^o| , pT] offer an interesting perspective 
for storing and manipulating arrays of atoms with the 
eventual possibility to scale, parallelize, and miniaturize 
the atom optics devices needed in quantum information 
processing. Moreover, optical microtraps can take ad- 
vantage of the fact that most of the current techniques 
used in atom optics and laser cooling are based on the 
optical manipulation of atoms . Many of the require- 
ments for the implementation of quantum computation 



p2| have been recently demonstrated in optical micro- 
traps containing ~ 100 atoms per site e.g., selec- 
tive addressing of single trap sites, and initializing and 
reading-out of quantum states in each site. In addition, 
the possibility to store and detect single atoms in optical 
dipole traps has been reported p7[ ]. 

With the demonstration of single-qubit gates being 
straight forward, what remains to be experimentally 
demonstrated is the capability of these optical microtraps 
to perform two-qubit quantum gates. The most promi- 
nent examples of such gates include the CNOT gate, the 
phase gate and the VSWAP gate ||, The latter 

transforms states |0)|1) and |1)|0), written in the com- 
putational basis, to maximally entangled states, while 
leaving |0)|0) and |1)|1) unaffected, in such a way that 



after the successive application of two vSWAP gates the 
states of the qubits are interchanged. Each one of these 
two-qubit gates, together with arbitrary single qubit op- 
erations, is universal, i.e., allows to perform any quantum 
algorithm. In practice, the particular two-qubit gate to 
be implemented will depend on the physical system under 
consideration. 

With respect to neutral atoms, several different phys- 
ical mechanisms to perform two-qubit gates have been 
proposed, ranging from cold controlled collisions |Q, ||] 
and dipole-dipole interactioiisj||, |l^, to purely ge- 
ometric quantum evolution p4[ . In the cold coUisional 
case, a two-qubit phase gate was proposed by adiabati- 
cally approaching two traps or by instantaneous state- 
selective switching of the trapping potentials § . In both 
cases the qubit was encoded in some internal degrees of 
freedom of the atoms, e.g., spin, Zeeman or hyperfine lev- 
els. For cold collisions to take place the atoms have to be 
brought to close distances, such that their quantum sta- 
tistical nature has to be taken into account. A detailed 
study of the role of the bosonic or fermionic character of 



2 



particles in the context of quantum information in atomic 
waveguide structures has been done by E. Andcrsson et 
al. Hi. 

Here we address the problem of implementing a quan- 
tum gate by adiabatically approaching two bosonic 
atoms, each stored in a different microtrap. In contrast 
to the proposals mentioned above, we assume the qubit 
to be implemented in the motional states of the atoms, 
i.e. an atom in the ground or the first vibrational state 
of the trap represents |0) or |1), respectively. Note that, 
as for the ion-trap case, the observation of neutral atoms 
cooled down to the ground and first vibrational states as 
well as superposition states in one-dimensional traps has 
been achieved ||2^. To perform the gate operation, we 
apply the steps outlined in Fig. 1. Initially, the two mi- 
crotraps are far apart such that the interaction between 
the two atoms is negligible. Then we adiabatically move 
both traps close together such that tunneling and cold 
controlled collisions become important. The dynamics of 
this process strongly depend on the particular motional 
state of the atoms and we can make use of this fact to 
control the interaction such that, after the eventual sepa- 
ration of the traps, the desired gate operation is realized 
with each trap again containing only one atom. 

To be more specific we consider here laser-cooled ru- 
bidium atoms stored in optical microtraps [ pT| , assuming 
that each trap contains initially only one atom. We will 
show that the V SWAP gate is the most natural quantum 
gate to be implemented when the qubit is encoded in the 
motional states of the atoms and interaction takes place 
through tunneling and cold collisions. This result applies 
to both, ®^Rb and ®^Rb, although they have negative and 
positive scattering length, respectively. In particular, we 
will demonstrate that a quantum gate of ~ 20 millisec- 
onds duration can be performed in state-of-the-art opti- 
cal microtraps. Very recently, E. Charron et al. have 
proposed the implementation of a phase gate in an opti- 
cal lattice with, as it is done also here, the qubit encoded 
in the motional states. In this case, a controlled interfer- 
ence set-up was proposed to perform a high-fidelity gate 
with operation time of 38 ms. 

The paper is organized as follows. In Section II, we 
introduce the physical model. Section III is devoted to 
the implementation of the VSWAP gate. In Section IV 
we discuss some practical considerations. And, finally, 
section V summarizes the results and presents the con- 
clusions. 



II. MODEL 

In this section we will first write down the Hamiltonian 
for the two atoms stored in the microtraps and discuss 
the interaction mechanism. We will introduce a time- 
dependent orthonormal set of single-particle states for 
each trap that is also orthogonal to the states of the other 
trap for arbitrary distances between the two traps. These 
single-particle states will make it possible to expand the 



wavefunction in a set of two-particle orthonormal states. 
This representation has two important advantages: (i) 
it allows to compute entanglement throughout the whole 
gate process; and (ii) it strongly reduces the computa- 
tional time required to simulate a quantum gate opera- 
tion with respect to a direct numerical integration of the 
Schrodinger equation for the two-particle spatial wave- 
function. Finally, we will discuss the physical implemen- 
tation of the qubits and its implications for the quantum 
gate operations. 



A. Hamiltonian 

The Hamiltonian governing the dynamics of the two 
atoms in a time-varying particle-independent trapping 
potential V(r, t) can be written as 
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+ U{ri- fa) , 



(1) 



where m is the mass of the atoms, f^ and pi are the 
(three-dimensional) position and momentum operators 
for atoms 1 and 2, and U {ri — fa) accounts for the in- 
teraction between the two atoms. 

To simplify the problem, we take the trapping potential 
shape to be time-independent along y and z directions: 



V{r, t) = v(x, t) + Vp{y) + Vp{z), 



(2) 



and assume much stronger confinement in y and z di- 
rections than in x, such that transverse excitations can 
be neglected. In fact, we will consider that both atoms 
are cooled down to the y and z vibrational ground states 
and remain there during all the interaction process. Ex- 
plicitly, we take the following one-dimensional potential 
to describe the two microtraps separated by a distance 
2a{t): 



v{x, t) 



{x + a{t)f 0{-x) + {x- a{t) f e{x) 



(3) 



where oj^ is the trapping frequency in the x direction, 
and 9{x) is the step function. 

The temporary variation of the trap distance is 
sketched in Fig. 1. Initially the traps are separated by a 
distance 2aj^g_^. The process of slowly approaching them 
to a minimum separation 2aniin takes a time tr and is 
modeled by the first quarter of a period of a cosine. Then 
we let the atoms interact for a time ti and, finally, we 
slowly separate the traps. 

For cold bosonic atoms, the dominant coUisional in- 
teraction is due to s-wave scattering which can be 
described by a contact potential of the form 
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FIG. 1: Separation of the traps as a function of time, tr 
and ti are the approaching/separating and interaction times, 
respectively. At Omax atoms located in different traps do not 
interact while at amin tunneling and cold collisions take place. 



where at is the s-wave scattering length of the rubidium 
atoms, e.g., in the spin triplet at — —369 aq for ^^Rb and 
at — 106 oo for ^^Rb with gq being the Bohr radius. As 
long as both atoms remain in the transverse vibrational 
ground states, we can integrate out the corresponding de- 
grees of freedom and obtain an effective one-dimensional 
interaction potential B 



u{xi — X2) — 2athujpS{xi — 



(5) 



where LUp is the transverse trapping frequency. Eqs. (j^) 
and (H) allow us to reduce the complexity of the problem 
to one dimension. 



B. Single- particle states 

We will implement the qubits into the ground and first 
excited vibrational states of each trap, i.e., we will use 
the motional states of the atoms. When the two traps are 
far apart, i.e., aa ^ 1 with 1/a = \/h/ muj^ being the 
position uncertainty of the ground state, these states are 
the energy eigenstates of two displaced one-dimensional 
harmonic oscillators: 



^1/4 



(6a) 

a(x±a), (6b) 



with L and R labeling the left and right trap, respectively. 
As we approach the two traps, these single-particle states 
overlap and are no longer orthogonal. To numerically 
integrate the Schrodinger equation and to compute en- 
tanglement throughout the gate process, we construct an 
orthonormal single-particle basis for arbitrary distances 
of the two traps by applying the Gram-Schmidt method 
(see Appendix A). If we denote these new single-particle 
states by with i = 0,1,2,3... and s = L,R then it 



holds s(*|i)t = SijSst- The four states that for large dis- 
tances correspond to the two lowest states of each trap 
read 
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where ^^(a) and S,t{a) are given in Eqs. (A5). For large 
separation of the traps, i.e. aa ^ 1, we have = ^~ 
for all i and thus the \i)L,R become the eigenstates of 
a single harmonic trap centered at =Fa. Notice that the 
|i)L._R, states have the following symmetry under parity 
transformation: {x\i)L,ji (~1)* (2;|i)_R.L- The general 
proof is given in Appendix A. This property obviously 
holds for the \i)L,R, and the \i)L.R. are constructed such 
that this symmetry is maintained. 

Although above we have written only four states, for 
all simulations using these orthogonalized states we will 
include all states up to \S)l,r- 



C. Two-particle states 

Let us motivate the two-particle basis which we will 
use. On the one hand, it must satisfy bosonic statistics, 
i.e. the basis states have to be symmetric under the per- 
mutation of the particles. On the other hand the Hamil- 
tonian of this system is symmetric with respect to parity 
transformation, i.e. H{x) = H{—x), and therefore does 
not couple states of opposite parity. For this reason we 
will introduce basis states with well-defined parity. If for 
this description we limit ourselves again to the four low- 
est single-particle states, then the bosonic two-particle 
sector forms a ten-dimensional Hilbert space. Here, we 
use the following notation |to(1))s 18) \n{2))t = \rn)s\n)t 
with I and 2 labeling the atoms and s,t — L,R. Thus, 
the bosonic two-particle basis reads: 



00) + 

01) + 

11) + 

00) + 

01) + 

n)+ 



;^(|o)Jo)^ + |o)^|o)J, 
^(|o)L|T)« + |T)^|o)i 

-|o)«|t)l-|t)l|o)J' 

;^(|T)L|T)« + |T)fi|T)J' 
;^(|o)l|o)l + |o)«|o)J' 

^(|0)l|T)l + |T)l|0)l 



V2 



{\1)l\1)l + \1)r\1) 



(8a) 

(8b) 
(8c) 

(8d) 

(8e) 
(8f) 
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and 



|oo)- 
|oT)- 

ITT)- 



+ |o)«|i)i + |i)Jo)J, 
;^(|o)l|o)l-|o)h|o)J 

^(|0)l|T)l + |T)l|0)l 
+ |0)JT)^ + |T)^|0)J. 

;^(|t)JT).-|t)h|tU 



(9a) 
(9b) 

(9c) 
(9d) 



The notation at the l.h.s. of Eqs. (||) and (H) means the 
foUowing: superscripts + or — indicate that the two- 
particle state has positive or negative parity, respectively, 
while the tilde accounts for states where, for aa 3> 1, 
both atoms are in the same trap, i.e., double-occupancy 
states. It is easy to check the symmetry of these two- 
particle states under the exchange of the atoms by mak- 
ing use of the parity property of states |^)^ ^ discussed 

after Eqs. (7). 

In addition it is worth to mention that in our simu- 
lations we will consider up to eight single-particle states 
which gives rise to a bosonic two-particle Hilbert space 
of 36 states (20 states having positive parity from which 
10 correspond to double occupancy; and 16 states hav- 
ing negative parity with 10 accounting for double occu- 
pancy). Finally note that the fact that we are able to 
expand the wave-function into this finite number of two- 
particle orthogonal states has also an important advan- 
tage with respect to the time needed for the simulation of 
a gate operation. We have checked the accuracy of the re- 
striction of the simulation to this subspace by comparing 
the results of the simulations to a direct numerical inte- 
gration of the Schrodinger equation for the two-particle 
spatial wavefunction which is about four orders of mag- 
nitude slower. 



D. Physical implementation 

We start from two well separated traps, each contain- 
ing one atom. In this situation we can neglect the bosonic 
nature of the particles and forget about the symmetriza- 
tion Only then it is possible to speak about well- 

defined qubits and we choose to introduce labels A and 
B for the two qubits by labeling the atom found in the 
left trap by A and the atom in the right trap by B. 

With the two traps far apart, single qubit operations, 
e.g., a Hadamard gate, can be realized by using two laser 
pulses in a Raman configuration focused solely on one 
of the traps. The quantum gate operation between two 
qubits is much more involved. As we approach the traps, 
due to tunneling there will be a non-vanishing probabil- 
ity to find both atoms in the same trap. Thus we can no 
longer distinguish the atoms such that bosonic statistics 



become important and the qubits are no longer well de- 
fined. If, however, we approach and separate the traps in 
such a way that finally there is again one atom in each 
of the well separated traps then we can attribute (new) 
labels A and B to them in the same way as before. 

These considerations suggest the following mapping of 
the states of the computational basis into the two-particle 
basis states of Eqs. (0) and (||): 
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- Ill) 



1 
1 

71 



(|01) + 
(|01) + 



101) 
101) 



(lOa) 
(10b) 

(10c) 
(lOd) 



Note that the two particle states at the r.h.s of ( pi]| ) have 
a trivial evolution at the trapping frequency (or multiples 
of it) that can be removed by including this phase in the 
definition of the single particle states. 

We will take states (10) as the starting set for the gate 
operation and, after setting the initial state, we will adi- 
abatically realize the gate. In this adiabatic regime, if 
we start in an energy eigenstate the system will follow 
this time-dependent energy eigenstate during the whole 
gate process. The only allowed transitions are those cor- 
responding to states that (i) are initially degenerate in 
energy, and, at short distances, (ii) become coupled via 
tunneling and/or cold collisions. Therefore, in order to 
find the most suitable gate to be implemented in this 
system, we have to identify these resonant couplings. 

For this aim we will first discuss the ideal case for which 
there is no interaction between the atoms, i.e. the case 
where at — in We then have the following resonant 
couplings: 



|00) + 
101) 

111)+ 



|00) + 
101)- 

111)+ 



1 10) - 1 10) 



101) 



(11a) 
(lib) 
(11c) 



where |01) = ^75(101)+ + |01)-) and |1G) = ^(|01)+ - 

|01)~). Therefore, there is a non-negligible probability 
(even if we move the two traps adiabatically) to have both 
atoms in the same trap after the gate operation. Note 
that the kinetic and trapping terms of the Hamiltonian 
do not directly couple |01) with |10) since they are single 
particle Hamiltonians and, therefore, they do not allow 
for the simultaneously change of the motional states of 
both atoms. The coupling between |01) and |10) is medi- 
ated through the double occupancy states |01) and |10). 
Clearly, in the non-interacting case, a quantum gate op- 
eration always has to face with double-occuppancy which 
makes the problem hard to handle. 

Figures 2(a) and (b) show, for a particular parameter 
set, the final state of the system after the whole process 
of approaching and separating the traps as a function of 
the scattering length. In (a) the initial state is |01) and in 



5 




20 40 60 80 100 120 140 



a/a. 



FIG. 2: Populations of the final state of the system af- 
ter adiabatically approaching and separating the traps as a 
function of the scattering length. The initial state is (a) 
|01) and (b) lll)"^, respectively. The parameter setting is: 



1.25 X 10* s" 



7.9 X 10^ 



1/a 



241 nm, 



5, aminO! = 1.99, Ulxtr = 70, aud LOxti 



(b) 1 11) . Although the scattering length has a constant 
value that depends on the atom under consideration, it 
is used in this plot as a free parameter to illustrate the 
double-occupancy problem. Notice that by changing LOp 
it is possible to tune the strength of the effective inter- 
action potential, cf. Eq. (||). Figs. 2(a) and (b) clearly 
show that, for at = 0, double-occupancy is indeed very 
important in the final state of the system. 

The problem of double-occupancy is naturally sup- 
pressed when one takes into account the interaction be- 
tween the atoms. In this case, double-occupancy states 
are no longer degenerate with single-occupancy states 
and we can neglect the probability to find double oc- 
cupation in the final state by adiabatically moving the 
traps. Thus, in the presence of interaction, the resonant 
couplings read 



101) 
111) + 



|10) 
|02) + 



(12a) 
(12b) 



(|0)l|2>«+|2);^,|0)^ + |0)^|2)^ + 



where_|02)+ - xr/L\"/R ' r/RT/L ' r/nr/L 

' L 1^^)^) /2- Notice that now the coUisional interaction 
term (0) allows for the simultaneous change of the mo- 
tional states of both atoms. The role of these couplings 
is clearly shown in Fig. 2. In Fig. 2(a), where the initial 



state is |01), double-occupancy populations in the final 
state start to decrease and eventually vanish as soon as 
the scattering length is increased. When the initial state 
is 111)"*", Fig. 2(b), double-occupancy also vanishes as 
the scattering length increases, but then the population 
of state |02)+ becomes important. 

Therefore, the coupling given in Eq. (12a) suggests the 
implementation of a VSWAP ga te, a s long as we are 
able to suppress/control coupling ( [l2b| ). The degeneracy 
between jll)"*" and |02) + can be broken, for instance, by 
taking an anharmonic trapping potential such that the 
vibrational frequencies are no longer equally spaced. In 
addition, it is possible to adjust the interaction time in 
such a way that, at the end of the gate operation, state 
|02)+ is not populated. In what follows, we will focus on 
this last possibility. 



III. VSWAP GATE 

The VSWAP gate has the following effect on the states 
of the computational basis: 



\Q)a\0)b 
|0)^|1)b 

|i>aI0)b 
|i>Ji)b 



\0)a\0)b 

^|0)^|l)s 
1 - i 



(13a) 
l)^|0)5(13b) 



|0)^|1) 



1 -i 
1 + z 



1 



|1U|1>B 



aI0)b(13c) 
(13d) 



It is straightforward to check that the successive appli- 
cation of two VSWAP gates exchanges the states of the 
qubits, i.e., C/swap = C^vswap ' ^ ^Jsw^- ^^^^'^ 
mentioned before, the VSWAP gate together with single 
qubit operations suffices to realize any quantum algo- 
rithm which is not the case for the SWAP gate it- 
self. A simple way to prove this, consists of showing that 
the universal CNOT gate can be obtained from VSWAP 
gates and single-qubit operations. In fact, a possible se- 
quence is (see Appendix B): 

C^CNOT = HAcrJ^^ciBU^gyfj^pa^U -^af (14) 

where Ha and aA,B are single-qubit operations. Addi- 
tionally, sequences involving single-qubit operation ex- 
clusively on one of the qubits, e.g., only on A, can be 
realized Hi. 



A. Gate simulation 

To simulate the gate operation, we have numerically 
integrated the time-dependent Schrodinger equation for 
the Hamiltonian given in (3) and (5) with the two- 
particle wave-function expanded in the previously intro- 
duced two-particle basis. Figures 3(a)-(c) show the re- 
sult of a VSWAP gate operation for a scattering length 
of at = 106 flo corresponding to ^^Rb atoms in the spin 
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FIG. 3: Simulated VSWAP gate operation for ^'^Rb with 
at = 106 ao- The rest of parameters as in Fig. 2. The initial 
state of the system is (a) |00)+, (b) |01), and (c) |11) + . We 
note that for ( b) the final relative phases of |01) and |10) are 
as in Eq. (|l3b|). 



FIG. 4: Simulated VSWAP gate operation for ^'^Rb and 
the following parameter values: uj^ = 1.25 x 10* s~^, ujp — 
1.6 X 10*^ s"\ 1/a = 244 nm, SmaxQ = 5, aminO: = 1.956, 
ujxtr = 77, ujxti = 97.2, and at = —369 ao. The initial state 
of the system is (a) |00)+, (b) |01), and (c) |11)+. 



triplet. The parameter setting is as in Fig. 2 and the 
initial state is (a) |00) + , (b) |01), and (c) The pa- 

rameter values are chosen to reproduce the gate operation 
given in (jl^) as well as to suppress the |02)+ population 
in the outgoing state of Fig. 3(c). Notice that states rep- 
resenting double occupation are populated at close dis- 
tances for all three cases. However, these populations 
vanish after the eventual separation of the traps since 
the traps are moved adiabatically and single-occupancy 
states are not degenerate with double-occupancy ones. 

For ^^Rb with negative scattering length it is slightly 
more involved to find parameters for the gate realization 
since, due to the attractive character of the interaction, 
double occupation states can more easily become reso- 
nant to single occupation states, e.g., |01)+ with |00) + . 
The parameters must be chosen to avoid this degener- 
acy between double and single occupation states. Fig. 4 
shows the result of a gate simulation for ^^Rb. Unlike 
for ^^Rb, Fig. 4(a) shows that starting from lOO)"*" state 
lOl)"*" is populated during the gate operation. 

On the other hand, it is important to notice that the 
results obtained for ®^Rb (Fig. 3) can be also directly 
implemented in ®^Rb by making use of the strong varia- 



tion of the scattering length in the vicinity of a magnetic 
field induced Feschbach resonance |Q . 

To check the accuracy of the previous simulations in 
which the two-particle wavefunction was expanded in a fi- 
nite set of states, we also have numerically integrated the 
Schrodinger equation for the two-particle spatial wave- 
function by using an operator split method and an FFT 
routine. Fig. 5 shows the results of this integration for 
the same parameter values as in Fig. 3 Q. The snap- 
shots give the joint-probability distributions for the two 
particles for three different initial states: (a) |00)+, (b) 
|01) and (c) 111)'*'. The bosonic nature of the atoms 
manifests in the symmetry of the joint-probability distri- 
bution along the diagonal xi — X2- In (a) and (c) the 
final state coincides with the initial one in accordance 
with Eqs. (p^). In (b) |01) evolves towards the maxi- 
mally entangled state [(1 -I- i)|01) + (1 - i)\lO)] /2 whose 
joint-probability distribution corresponds to the donut- 
like shape of the last frame. 

The accuracy of the simulated gate operation U with 
respect to the perfect gate operation U y' g^Ap as given 
by Eqs. ( p^ ) is computed through the averaged fidelity. 
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FIG. 6: Averaged fidelity of tlie gate operation in tlie parame- 
ter plane tr versus amin- The interaction time is (a) uixti — 69 
and (b) tUxti = 20. The rest of the parameters as in Fig. 3. 



i.e. 



F = Tr[[/pC/t PC/^SWAP ] (15) 

where the average is taken over the four orthogonal pure 
input states p from Eqs. (jl^). Figure 6(a) shows for 
®^Rb the averaged fidelity F of the gate in the parameter 
plane versus amin- The rest of the parameters are as 
in Fig. 3. Clearly, the fidelity is very sensitive to the 
minimum distance due to the exponential dependence of 
tunneling at this distance. Note that the fidelity of the 
gate operation corresponding to the parameters of Fig. 3 
is > 0.9997 with a gate duration of 2tr + U ^ 17 ms 
for Lj^ = 1.25 X 10^ s-\ 

An important issue is how much the gate duration 
can be decreased while maintaining a high fidelity. In 
Fig. 6(b) the gate duration is reduced by a factor of 
two which reduces the fidelity by around one order of 
magnitude. In fact, as soon as the rising time t^. is de- 
creased, non-adiabatic effects occur which in turn result 
in the population of several unwanted states, e.g. double- 
occupation, in the final state of the system. However, it 
could be possible to use the techniques developed in |3|] 
to optimize the speed of the gate operation, while sup- 
pressing excitations to these unwanted states. 



FIG. 5: Snapshots of the spatial two-particle wavefunction 
\'ip{xi,X2)\'^ for ^'^Rb. The parameters are as in Fig. 3. The 
horizontal and vertical axes of each plot show the coordinate 
of the first and second particle, respectively. Initially there 
is one particle in each of the traps. The symmetry along the 
diagonal xi = X2 is due to the bosonic statistics. A particle in 
the ground state of one of the traps corresponds to a gaussian 
distribution in the direction of the respective axis while one 
node corresponds to the first excited state. Thus the initial 
states are (a) |00) + , (b) |01), and (c) |11) + . The time for the 
snapshots is shown in (d). See for animated illustrations 
of the gate operation. 



B. Quantum correlations 

Let us briefly consider entanglement in the context of 
the VSWAP gate. As already discussed the initial and, 
as long as double occupancy is suppressed, also the final 
state consist of well separated and thus for practical pur- 
poses distinguishable particles, such that the usual no- 
tions of entanglement apply. Let us consider in particular 
the separable state |01) that under the VSWAP opera- 
tion evolves to the maximally entangled state from Eq. 
(131:). During the operation bosonic statistics become 
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important and thus we have to distinguish between sta- 
tistical correlations arising from symmetrization on one 
hand and quantum correlations which are useful in the 
context of quantum information on the other hand. Such 
a distinction has been discussed in j^, |3^ for fcrmionic 
two-particle states where the notion of Slater rank was 
introduced and a fermionic measure of entanglement was 
derived. In Q these methods have also been used to 
study correlations in the context of a quantum gate op- 
eration for two electrons in quantum dots. They have 
been translated to bosons in p4| and moreover a bosonic 
von Neumann entropy has been defined in | |35[ |. Here we 
want to study to which extent these techniques can be 
applied to our particular gate operation. 

Let us write a general two-boson state in an 
iV-dimensional single particle space as \v) = 
'^fj^iVijblbjl^l) where 6| and bi are bosonic cre- 
ation and annihilation operators and \ is the vacuum 
state such that bl\n) = The complex symmetric 
matrix Vij — vji is normalized as t'!:{v^v) — 1/2. If 
new bosonic annihilation operators b'i — Uijbj 
are introduced by a unitary transformation U of the 
single particle space, then v transforms as UvU'^ . Now 
we find that for every symmetric complex matrix v 
there exists a unitary U such that UvU^ is diagonal, 
i.e. UvU'^ = diag[Ai, . . . , Ar, 0, . . . , 0] with > 
p^ , |35[ |. Here r is called the Slater rank of \v) and 
\v) — Xib'lb'l\Q) its Slater decomposition. 

In our case the initial state |01) has Slater rank two 
while the final state U ^ gwAp l*-*-*-) Slater rank four. 
A bosonic measure of entanglement has been defined in 
Q only for the case of a two-dimensional single parti- 
cles space, but in a bosonic version of the von Neu- 
mann entropy has been proposed that, as a function of 
the Slater coefficients A,- reads 



(16) 



Sb ranges from Sb — for states with Slater rank 1 
to Sb = log2(-/V) for states with Slater rank N and 
all \i equals. For the case considered here we have 
Sb{\01)) = 1 and SB{U^g^^\01)) = 2. The function 
SB{t)/2 is plotted in Fig. 6 for the complete VSWAP 
operation (the factor is such that Sb for N = A ranges 
between and 1). The figure also shows the product 
'S' • Psingic where S is the von-Neumann entropy calcu- 
lated by projecting the state onto the space spanned by 
the set {|00)+, |01)+, |01)-, |11)+} and renormalizing. 
If now states \i)L are considered as being distinguishable 
from states \i)R then S can be calculated as for distin- 
guishable particles. Psingic is the probability to find the 
state in the space spanned by the given set. 

Fig. 7 shows that in the limit of large spatial separation 
of the particles neither the Slater rank nor the bosonic 
von Neumann entropy from p5| coincide with the notion 
of Schmidt rank and von Neumann entropy for distin- 
guishable particles. There are two reasons for this obvi- 
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FIG. 7: Correlations for the VSWAP gate operation applied 
to |01). Bosonic von Neumann entropy Ss/i (solid line) and 
von Neumann entropy S (dashed line) after projection onto 
the space spanned by {|00)+, |01)"'', jOl)", |11) + } multiplied 
by the probability psingic to find the state in the space spanned 
by this set. 



ous difference: (i) Slater rank and bosonic von Neumann 
entropy are invariants under arbitrary transformations of 
the single particle space, while in our case the partition 
of the full single-particle space into the space spanned by 
\i) L on the one hand and \i)ii on the other is rather fixed; 
(ii) they do not take into account that here we want to 
avoid double-occupation in the final state. 

In P. Zanardi follows a different approach to cal- 
culate an entropy for indistinguishable particles which 
consists of mapping the bosonic (or fermionic) state to 
an occupation number basis and calculating the von Neu- 
mann entropy in this basis. If this is applied here then 
the entropy raises from for |01) to 1 for U ^ swAp l^-*-) 
but in between it has a maximum of ~ 1.4. As can be 
seen from the analysis in this is due to the presence 
of spatial entanglement during the operation. It is how- 
ever not clear whether this type of entanglement is useful 
in this context where the two qubits are explicitely im- 
plemented in the vibrational states of each trap. These 
results thus call for further investigation of the problem 
of quantum correlations in bosonic states. 



IV. PRACTICAL CONSIDERATIONS 

Scalable systems of optical microtraps based on the 
dipole force can be realized by focusing a single red- 
detuned laser beam with a microlens array |p^ , pl]| . The 
temporal evolution of the trap separation as shown in 
Fig. 1 can be realized (i) by using two parallel laser beams 
focused in such a way that the trapping potentials are 
longitudinally shifted along the common laser beam di- 
rection or (ii) by illuminating the microlens array with 
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two laser beams under slightly different angles with the 
possible inclusion of an additional two-dimensional con- 
fining potential perpendicular to the direction of the trap 
displacement [ |l9| . For the parameters we used in the 
gate simulations, the minimum distance of aamm ^ 2 
corresponds to a separation of the traps of 1 um which 
is achievable in the present optical microtraps pl| . With 
laser powers of 1-10 mW per trap, rubidium atoms can be 
trapped with typical trapping frequencies along the laser 



beam direction of < 



10 -10 s while the transverse 



trapping frequencies can be one or two orders of mag- 
nitude larger |l^. Additionally, sideband cooling could 
be applied to cool the atoms to the ground state of each 
trap in all dimensions. 

In the optical microtrap experiments, the trapping po- 
tential is gaussian-shaped with typical depths of 1-10 
mK X kB jl^. For a single trap it is thus a good ap- 
proximation to assume a harmonic potential for the lower 
lying states. For two traps being close together the ac- 
tual potentials will deviate from the form assumed in (|^). 
Nevertheless, it is possible to generalize the methods ap- 
plied here to these particular potentials. To allow for 
single site addressing, typical initial distances between 
the traps have to be about 5-20 /im which for rubidium 
means aamax — 10 — 40 instead of aamax = 5 which 
we used in our previous calculations. In this case, it is 
straightforward to estimate that the process of adiabat- 
ically approaching the traps to the interaction distance 
and eventually separating them again requires another 
~ 10 ms. Therefore, the whole gate process could be 
realized in a few tens of milliseconds which is enough 
for a proof-of-principle experiment. In fact, the lifetime 
of the atoms in the traps is about 100 — 1000 ms while 
spontaneous scattering occurs in ~ 10 ms. However for 
strongly confining trapping potentials only a small frac- 
tion of the occurring spontaneous scattering processes 
leads to a change in the vibrational state |Q . 

In a recent paper, E. Charron et al. [ p8| proposed the 
realization of a phase gate in an optical lattice where 
the qubits were also implemented in the motional states. 
Two linear counterpropagating beams from the funda- 
mental and first harmonic of a CO2 laser were used to 
produce an intensity gradient optical lattice. The barrier 
between two neighboring traps could be raised or lowered 
by changing the intensity ratio between the two beams. 
We notice that the realization of a vSWAP gate as dis- 
cussed here should also be possible in this setup although 
the implementation in optical microtraps presents some 
advantages such as being not sensitive to the phase fluc- 
tuations of the lasers. 



V. CONCLUSION 

We have investigated quantum computation in optical 
microtraps with the qubits implemented in the motional 
states of neutral atoms, and tunneling and cold controlled 
collisions accounting for the interaction between two dif- 



ferent qubits. A time-dependent two-particle orthogonal 
basis has been introduced to simulate the gate opera- 
tion and to compute entanglement throughout the whole 
gate process. The bosonic statistic nature of the parti- 
cles and its role in entanglement has been discussed in 
detail. We have demonstrated the capability of optical 
microtraps to realize a high-fidelity ySWAP gate opera- 
tion in the few tens of milliseconds range. Finally, some 
practical considerations for the physical implementation 
of this quantum gate have been discussed. 
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APPENDIX A: GRAM-SCHMIDT 
ORTHONORMALIZATION 

In this appendix we will show how to construct the 
time-dependent orthonormal single-particle states, de- 
noted by 1?)^ and from the harmonic oscillator en- 
ergy eigenstates and for the left and the right 
trap, respectively. We start by defining states involving 
one state of each trap: 

1 



^7iLN)L±(-iri 



i = 0,1,2,3,, 



(Al) 



where the superscript + (— ) indicates positive (nega- 
tive) parity with respect to the middle between the two 
traps. We then group these states according to their 

parity in two sets = ||0)* , |1)* , ...| and focus first 

on the positive parity set . This set contains states 
that are neither orthogonal nor normalized. To perform 
the orthonormalization, we use the Gram-Schmidt (GM) 
method starting with the following normalized function: 

{x |0) + 



/ I {x |0)+ \ dx 



(A2) 



Then, we define the first linearly independent function 



as: 



cjy[{x,a) 



(x ID 



-aio0(|(a;, a) 



/|(a;|l) + aio(l)Q{x,a)\ dx 



(A3) 
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where aio = — / a)(a; |l)~''da; which guaranties 

{4>o \4>i) = 0. We repeat this procedure to obtain 
the rest of the hnearly independent functions (j>2,<j)^, ■■■ 
with positive parity. In an analogous way, we deter- 
mine from S~ the set of hnearly independent functions 
jcf)^ ,cj)2 ,(j>^ , An important feature of the GM 
method when applied to a set of states with the same par- 
ity is that the constructed orthonormal states retain the 
parity of the original set of states. Thus states from } 
and {(pf} have positive and negative parity, respectively 
and, therefore, the whole set {^q , </)f , (/)^, 0^, ...} is or- 
thonormal. Explicitly, the first four orthonormalized 
functions read 



v3/2 



where 



1 



(A4a) 
{x |0)fi^W 



(A5a) 



ga^a^ ± 1) fe-'"' - e-"^-' 



4A5b) 



For the sake of brevity, we do not explicitly show the 
analytical expressions for the rest of the 0^ . 

Once we have obtained the orthonormal set {ipf}, it 
is straightforward to write down the single-particle basis 
that we will use: 



{x 



(A6a) 
(A6b) 



These states are orthonormal due to the orthonormality 
of the (pf and in the limit aa 3> 1 become the corre- 
sponding harmonic oscillator energy eigenstates for each 
trap. These new orthonormal states do not have in gen- 
eral a well defined parity with respect to the center of 
the corresp ond ing trap but it is straightforward to check 
from Eqs. ( |A6D that they satisfy the following property 
under parity transformation with respect to the middle 
of the traps: 



{x\i)L,R (-1)' {x\i)Rj 



(A7) 



APPENDIX B: UNIVERSALITY OF THE VSWAP 
GATE 

Our goal here is to write down the sequence of steps 
required to build the CNOT gate, which, in the compu- 



tational basis \0)a\0)b, \0)a\1)b, \1)a\Q)b and \1)a\1)b, 
reads 



CNOT — 



10 
10 
1 
10 



(Bl) 



from the VSWAP gate 



U. 



VSWAP 



/ 1 \ 

^ ^ 

1 



\0 



(B2) 



The single qubit operations we need are on one hand the 
Hadamard gate 



H 



V2 li -ly 



(B3) 



and on the other hand the following combination of iden- 
tity and Pauli CTj matrices: 



1 

-i 



(B4) 



Let us call Ha,b and aA,B the corresponding single-qubit 
operations for qubit A or B. Now it is easy to check that 
the following combination of single-qubit operations and 
VSWAP gates yields the phase gate: 



PHASE 



,-1 



(B5) 



/ 1 
10 
10 

\ -1 



This sequence is not unique, and more sophisticated se- 
quences involving single-qubit operations on only one of 
the qubits can be implemented . Finally, to obtain the 
CNOT gate it is enough to apply a Hadamard gate on 
the qubit A at both sides of Eq. 



U, 



CNOT 



— HaUphaseHa- 



(B6) 
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